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spheres 

Andres Santos, 1 '*) Gilberto M. Kremer 2 '**) and Vicente Garzo 1 '***) 

1 Departamento de Fisica, Universidad de Extremadura, E-06071 Badajoz, Spain 
2 Departamento de Fisica, Universidade Federal do Parana, Curitiba, Brazil 

The aim of this work is to explore the combined effect of polydispersity and roughness 
on the partial energy production rates and on the total cooling rate of a granular fluid 
mixture. We consider a mixture of inelastic rough hard spheres of different number densities, 
masses, diameters, moments of inertia, and mutual coefficients of normal and tangential 
restitution. Starting from the first equation of the BBGKY hierarchy, the collisional energy 
production rates associated with the translational and rotational temperatures (T* r and T/ ot ) 
are expressed in terms of two-body average values. Next, those average values are estimated 
by assuming a velocity distribution function based on maximum-entropy arguments, allowing 
us to express the energy production rates and the total cooling rate in terms of the partial 
temperatures and the parameters of the mixture. Finally, the results are applied to the 
homogeneous cooling state of a binary mixture and the influence of inelasticity and roughness 
on the temperature ratios T{ T /T{ ot , Tj'/T* 1 , and T 2 rot /T[ ot is analyzed. 

§1. Introduction 

A granular fluid is usually modeled as a system of identical, inelastic smooth hard 
spheres with a constant coefficient of normal restitution a. Despite its simplicity, 
this model has been useful to capture the basic properties of granular flows. 1 - 1 On the 
other hand, the model can be made closer to reality by introducing more ingredients, 
such as coefficients of normal restitution depending on the impact velocity, 2 ) presence 
of an interstitial fluid, 3 ) non-spherical shapes, 4 ) polydispersity, 5 ) or roughness. 6 ) Of 
course, the few citations in the preceding sentence are just representative of many 
works reporting features not accounted for by the simple monodisperse smooth- 
sphere model. 

In this paper we will focus on the two latter ingredients, namely polydispersity 
and roughness. These properties are especially relevant, not only because beads and 
grains are unavoidably polydisperse and rough, but also because any of these ingre- 
dients unveils an inherent breakdown of energy equipartition in granular fluids, even 
in homogeneous and isotropic states. In the case of multi-component granular fluids, 
most of the studies have considered the inelastic smooth-sphere description. Some 
of the problems addressed include non-equipartition in homogeneous states, 7 )~ 10 ) 
Navier-Stokes transport coefficients, n )~ 15 ) and segregation phenomena. 16 )~ 20 ) 

Concerning the case of inelastic rough spheres, most of the works we are aware 
of are restricted to monodisperse systems. 21 )~ 41 ) Analogously to what happens with 
the coefficient of normal restitution a, the simplest model accounting for friction 
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during collisions assumes a constant coefficient of tangential restitution (3. While 
a is a positive quantity smaller than or equal to 1 (the value a = 1 corresponding 
to elastic spheres), the parameter (3 lies in the range between —1 (perfectly smooth 
spheres) to 1 (perfectly rough spheres). Except for a = 1 and /3 = ±1, the to- 
tal kinetic energy is not conserved in a collision. Some of the early attempts to 
develop a kinetic theory for rough spheres were carried out by Jenkins and Rich- 
man 21 ) and Lun and Savage, 22 ) who applied their approaches to the simple shear 
flow problem. The influence of roughness in shear flows has also been studied by sev- 
eral authors, 23 ^ -29 ) usually assuming that the spheres are nearly smooth and nearly 
elastic. In an extensive paper, 30 ) Goldshtein and Shapiro obtained the collisional 
energy production rates associated with the translational and rotational degrees of 
freedom by using Maxwellian forms for the distribution functions. The result was 
applied to the evaluation of the ratio between the translational (T tr ) and rotational 
(T rot ) temperatures in the homogeneous cooling state (HCS). The time evolution 
of the ratio T tr /T rot towards its HCS asymptotic value has been widely analyzed, 
both theoretically and by means of molecular dynamics, by Luding, Zippelius, and 
co-workers. 31 ) -37 ) Other studies involving roughness include vibration with rough 
walls, 38 ) a micropolar fluid model for granular flows on a slope, 39 ) derivation of hy- 
drodynamic constitutive equations from the Boltzmann equation for nearly smooth, 
nearly elastic granular gases, 40 ) and correlations between the rotation axis and the 
translational direction. 41 ) 

The studies about multi-component rough-sphere systems are much scarcer. To 
the best of our knowledge, only the case of a fixed particle immersed in a bath 
of thermalized point particles has been addressed. 42 ) -44 ) On the other hand, the 
general case of a polydisperse system made of mobile particles of different coefficients 
of normal and tangential restitution (q^- and fiij) has not been studied yet. In this 
paper, we address one of the basic aspects of the problem, namely those related to the 
partition of the total energy. In order to characterize the effect of collisions on energy 
partition, we focus on the partial productions rates £*j and measuring the rate of 
change of the translational and rotational kinetic energies, respectively, of particles 
of component i due to collisions with particles of component j. A combination of 
and Qj l gives the total cooling rate ( of the mixture. 

Starting from the collisional rules worked out in Section 2, the collisional rates 
of change of momentum (linear and angular) and energy (translational and rota- 
tional) are expressed as linear combinations of two-body average values in Section 
3. These averages are evaluated in terms of the partial temperatures T* r , T/ ot , T* r , 
and Tj ot by assuming a maximum-entropy two-body distribution in Section 4. The 
expressions for and ( obtained in Section 5 extend previous results derived 

for monodisperse rough spheres 30 )' 33 ) and for polydisperse smooth spheres. 7 ) An 
application to the HCS of a binary mixture is carried out in Section 6 with some 
illustrative examples. Finally, some concluding remarks are presented in Section 7. 
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§2. Collisional rules 

Let us consider the collision between two hard spheres of masses and rrij, 
diameters U{ and o~j , and moments of inertia Ii and Ij . The latter two quantities can 
be equivalently characterized by the dimensionless parameters 

-- --- J (2-1) 



9 5 ,IJ 7 — o * 



The value of Ki depends on the mass distribution within the sphere and runs from 
the extreme values Ki = (mass concentrated on the center) to Ki = | (mass 
concentrated on the surface). If the mass is uniformly distributed, then Ki = |. Let 
us denote by Vjj = Vj — Vj the pre-collisional relative velocity of the center of mass 
of both spheres and by and Uj the respective pre-collisional angular velocities. 
This is sketched in Fig. 1, where a = (ta — ) / 1 r ^ — | is the unit vector pointing 
from the center of % to the center of j. The velocities of the points of the spheres 
which are in contact during the collision are 

Wj=Vj - yffXWi, Wj = Vj + yO- X Ufj, (2-2) 



the corresponding relative velocity being 



wy = Vij -ax Sij, Sij = ycjj + ( 2 ' 3 ) 



Post-collisional velocities will be denoted by primes. Conservation of linear and 
angular momenta yields 37 ) 



m 



v. + rrijv'j = rriiVi + rrijVj, (2-4) 

<Vj. 

2-5) 



JjUJ- -mi^xvj = IiUi -rrii^S xv„ IjJj + mj^-ax v'j = IjU)j + mj^-ax Vj . 
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Equations (2-4) and (2-5) imply that 



v- = Vj-^-Qy, Vj = Vj + — Q i: j, (2-6) 



m,; 



where Qij is the impulse exerted by particle i on particle j. Therefore, 

y'ij = Vij - ^-Qij, w'ij = Wij - -^-Qij + — i— a x (a x Q^) , (2-8) 



where 



nii + mj Kinii + Kjmj 



(2-9) 



are the reduced mass and a sort of reduced inertia-moment parameter, respectively. 

To close the collisional rules, we need to express Qij in terms of the pre-collisional 
velocities and the unit vector <r. To that end, we relate the normal (i.e., parallel to 
<x) and tangential (i.e., orthogonal to a) components of the relative velocities 
and w^- by 

a ■ w'ij = -oiijd- ■ Wij, & x w'ij = -Pijcr x Wy. (2T0) 

Here, aij and fyj are the coefficients of normal and tangential restitution, respec- 
tively. The former coefficient ranges from = (perfectly inelastic particles) 
to = 1 (perfectly elastic particles), while the latter runs from fiij = — 1 (per- 
fectly smooth particles) to fiij = 1 (perfectly rough particles). Inserting the sec- 
ond equality of Eq. (2-8) into Eq. (2T0) one simply gets <x • Qij = a^a • w^- and 
a x Qij = fiijB x Wij, where 

aij = mij(l + aij), ^ = ^1(1 + ^). (2-11) 

J- ~t~ nij 

Therefore, 

Qij = oiij{vij ■ a) a + [vy -ax S^- - (v„ • a)S] , (2-12) 

where use has been made of the mathematical property <x x (<x x A) = (er ■ A)<r — A. 
Note that in the special case of perfectly smooth spheres {f3ij = — 1) one has f3ij = 0, 
so that er x Qij = 0. In that case, according to Eq. (2-7), the angular velocities of 
the two colliding spheres are unaltered by the collision. 
The total kinetic energy before collision is 

E i3 = ^vf + ^v] + \^ + kjj. (2-13) 
It can be checked (see the Appendix) that 

E'ij - Eij = (i _ f ) [{a x Vo .)2 + {9 x 8^)2 + 2{B x Vij) . Sij] 

-^(i-4)(--^) 2 - ( 2 - 14 ) 
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The right-hand side is a negative definite quantity. Thus, we observe that energy 
is conserved only if the particles are elastic (aij = 1) and either perfectly smooth 
(Pij = —1) or perfectly rough (fy = 1). Otherwise, E[- < Ey and kinetic energy is 
dissipated upon collisions. 

Equations (2-6), (2-7), and (2-12) give the direct collisional rules. The restituting 
collisional rules are 



v " = Vi - —Qij, v^ = v, + — Q: 



nii 



U ' " 2T ? X Qi ^' J i =Ui 3- ^T d x Qm 



K)" 



2U 



• j ' 



where 



a «3 ~\-~ | Pij 



Qij = ^(v.j ■ $)S + [Vij -ax Sij - (Vy • a)a] . 

Cxi - 



(2-15) 
(2-16) 

(2-17) 



Hj Pij 

Here the double primes denote pre-collisional quantities giving rise to unprimed quan- 
tities as post-collisional values. The modulus of the Jacobian of the transformation 
between pre- and post-collisional velocities is 



d(vi,u)i,Vj,u>j 



(2-18) 



§3. Collisional rates of change 



Let fj?\ri,Vi,u>i;rj,Vj,u)j;t) be the two-body distribution function with the 
normalization condition 

J dr., j dvi J dwi j drj J dvj J dcj j f^\r i ,v i ,u i ;r j ,v j ,u j ;t)=N i N j , (3-1) 

Ni being the number of spheres of component i. The one-body distribution function 
is 

fi(ri,Vi,LJi;t) = Nj 1 j drj J dvj J du>j f\f(r h v it uf, r jt v j: uj; t). (3-2) 

The marginal distribution functions associated with the translational and rotational 
degrees of freedom are 

fl x {r i ,\ i ;t) = J dWifi(ri,Vi,Ui;t), f- ot (r; , Ui ; t) = J dv, /jfo, v,, a;,; t). (3-3) 

Given a one-body function ^;(vj,u>j), we define its average as 

(Y>(V;,U>;)) = ^- J dVi J dWi1p(Vi,U>i)fi(Vi,U>i), n i = j dv i J dUifi(vi,Ui), 

(3-4) 
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where n{ is the number density of component i and, for the sake of brevity, we have 
omitted the spatial and temporal arguments. In particular, one can define partial 
temperatures associated with the translational and rotational degrees of freedom as 



where 



u 



Ei rrnni(\i) 



(3-5) 



(3-6) 



Ei m i n i 

is the flow velocity. Note that in the definition of JJ ot we have not referred the 
angular velocities to any average value because of the lack of invariance under the 
addition of a common vector to every angular velocity. The global temperature is 

i 

where n = J2i n i is the total number density. 

By starting from the Liouville equation and following standard steps, one can 
derive the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. 45 ) The 
first equation of the hierarchy reads 

dtfi(ri,Vi,u>i;t)+\i • V/ifc, v*, Ujji) = ^ J^, v^o^l/^], 

3 

where 

Jij[ru Vj, uf, t\flf] = afj J dvj J du>j J da 6>(v^- • S)(v tJ • a) 



(3-f 



1 



(2) 



(3-9) 



is the collision operator. Here, u^- = (crj + (Tj)/2 and use has been made of Eqs. 
(2-10) and (2-18). 

Multiplying both sides of Eq. (3-8) by ^(vj,cjj) and integrating over Vj and u>j 
one gets 



dtni(ip(vi,u)i)} + V • ni{viip(-Vi,u;i)} = ^2Jij[ip{vi,^i)\f, 



(2)i 
ij J 



(3-10) 



where 



^■[^(vi^i)!/^] = J dvi J dwi^{yi,u3i)Jij\yi,uii\flf] 

= afj J dvi J du>i j dvj J duj J da 0(v^ • a)(vij ■ a) 



flf{Ti,\i,uji;Ti + Gija,\j,ujj) [^(v-,w-) -^(Vi^i)] . 



(3-11) 
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Thus, nJ L Jij[i[)(vi,u>i)\flj ] is the rate of change of the quantity ip(vi,u)i) due to 
collisions with particles of component j. This rate of change is a functional of 

(2) 

the two-body distribution function f-j, as indicated by the notation, and it is in 
general a rather intricate quantity. The physically important cases are ^(v^ojj) = 
{mjVj, IiU)i, rriivf , hoof} . The corresponding rates of change are obtained by inserting 
Eqs. (A-5)— (A-8) into Eq. (3T1). Note that so far all the results are formally exact. 
To proceed, let us make the approximation 



(Vi,Wi)l4 ?> ] 



(3-12) 



where 



J da &(vij ■ ?)(vij ■ a) 



(3-13) 



j au <oyxij ■ (J)\\ij ■ 

(2) 

is the orientational average of the pre-collisional distribution f-j ' . Thus, Eq. (3T2) 

(2) 

replaces a detailed functional of f-j by a simpler one where the solid angle integral 
dveivij • ?)( Vij • a) [^(vJ.wO -^(vi,Wi)] (3-14) 



(2) 

can be evaluated independently of f-- . It is important to borne in mind that the 



approximation (3T2) is much weaker than the approximation fffl ~ f\t. } ■ Notwith- 



K2) 



standing this, the equality f\, = f^ ! holds if (a) the system is homogeneous and 



F(2) 



isotropic (regardless of the reduced densities n^af and njO~j), in which case / 



f and njO~j), in which case fjp 
only depends on |rj — r^- 1 , or (b) the system is in the Boltzmann limit (n^af — > 0, 

n j a j ~ y 0)' ™ which case one can formally take o~ij — > in the contact value of f^\ 
Therefore, the approximation (3T2) is justified if the density of the granular gas 
and/or its inhomogeneities are small enough so the value of fjp at contact is hardly 
dependent on the relative orientation of the two colliding spheres. 

In the remainder of this Section we particularize to ip(vi, 0Ji) = {mjVj, ijWj, mivf,Iii 
and express the rates of change n~ x Jij\}\){y Wj)] in terms of two-body averages of 
the form 

(^(vi,o;j; Vj^j)) = J dwi J dui J dvj J dk)j A(\i,u)i;Vj,u)j) 



(3-15) 



The results are (see the Appendix) 



71 i 1 Jij[ m i v i] 



2 / OLij + fyj 



(vy x Sij) 



(3-16) 



1 Jij[IiUi 



2 n a 
-njOijVi-P, 



3(%-Sy) - (v^ (Vij ■ S, ; )v,,) 



(3-17) 



A. Santos, G. M. Kremer, and V. Garzo 



Table I. Estimates of the two-body averages appearing in Eqs. (3-16)-(3-19), as obtained from the 
replacement (4-2). 



Quantity 



Estimate 



(VijVij) 

{vij x Sij) 
(vijSij) 

("y^.J ■ Si 3 -)Vij) 
(VijVi ■ Vij) 

(Sij ■ (v, xvj)) 
(VijSij) 

(vy ■ S«) 2 > 
(vijUi ■ Sij) 
(vT. 1 (vij ■ Sij) (v^ • Wi)> 

(Vij) 

(vij) 



5 (Oifli + (Tjflj) (Vij) 

i (aifli+ajflj)^) 

r tr 



TTITT TV 





(v 



131 



6Tf 



+ • 17j ) (v^) 

1/2 



2 \/ ( „L + J. 



i x . . ( z£ + 21 

■K \ TTli 7TLA 



3/2 



(aij + Aj) (VijVi • Vy) + ^"(Sij • (Vj X Vj)) 



~2 , o2 oo2 o2 



,(3-18) 



n- ^[i^ 2 ] = -njaf^ijl 3o-i(vijUi ■ Sy) - cr^^ 1 (v y - • Sy) (Vy • 



TTliKi 



2{v%) + Z(vijSl)-(v-Mwij-SijY) 



(3-19) 



§4. Estimates of the average values 



Equations (3-16)— (3-19) express the collisional rates of change of the main quan- 
tities as linear combinations of two-body averages of the form (3-15). They are local 
functions of space and time and functionals of the orientation-averaged pre-collisional 
distribution While, thanks to the approximation (3-12), Eqs. (3-16)-(3-19) are 
much more explicit than the exact results obtained from Eq. (3-11), they still require 
the full knowledge of f-j\ 

Suppose, for simplicity, that (vj) = (vj) = u and define the average angular 
velocities 

(u>i) = fii, (u>j) = flj. (4-1) 

-(2) 

Now, let us imagine that, instead of the full knowledge of , we only know the local 
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values of the two densities (rij and rij), the two average angular velocities ( fli and 
f2j), and the four partial temperatures (T* r , T- rot , T* r , and T* ot ). The question we 
want to address in this section is, can we get reasonable estimates of the two-body 
averages appearing in Eqs. (3-16)-(3T9) by expressing them in terms of rii, rij, f2i, 
f2j, Tj tr , T/ ot , Tj r , and TJ *? In the absence of further information, the least biased 
estimates are obtained from the replacement 



fij (Vi, WiJ V,- , Wj ) -> Xij I 47r2T tr^tr I eX P "^^^t? m 3 ^ 

xfi*(*H)f?*(!*j), (4-2) 

where Xij 1S the contact value of the pair correlation function. Equation (4-2) can 
be justified by maximum-entropy arguments, except that here we do not need to 
assume a Maxwellian form for the rotational distributions, given that the angular 
velocities only appear linearly or quadratically in Eqs. (3-16)-(3-19). Like in the 
approximation (3-12), it is important to stress that we are not making the strong 

~(2) 

claim that f>- (vj, u>f, Vj,u)j) is well approximated by the right-hand side of Eq. (4-2) 
[see Ref. 41)] but only the wekaer one that the two-body averages can be estimated 
by performing such a replacement. Table I displays those estimates. 

Although (vj) = (vj) has been assumed in the results shown in Table I, the 
generalization to (vj) / (vj) can be carried out following the same steps as done in 
Ref. 46) for smooth spheres. 

§5. Energy production rates and cooling rate 

The most characteristic feature of a granular gas is the energy dissipation taking 
place after each collision. In the model of inelastic rough hard spheres this is clearly 
apparent from Eq. (2-14). On the other hand, any of the four partial kinetic energies 
contributing to in Eq. (2-13) can either increase or decrease after a given collision. 
To characterize this effect at a statistical level, it is convenient to introduce the rates 
of change of the partial temperatures T^ T and Tf ot due to collisions of particles of 
component i with particles of component j. More explicitly, we define the (partial) 
energy production rates £*j and as 

^S^-^f^Wvi-u) 2 ], ^^-^^JiAh^Y (5-1) 

When collisions of particles of component i with all the components are considered, 
we get the (total) energy production rates 

Finally, the net cooling rate is 

^-K§) con =^ iT ™' +m - (5 - 3) 



,2 
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In contrast to the energy production rates denned in Eqs. (5-1) and (5-2), the cooling 
rate £ is positive definite, i.e., collisions produce a decrease of the total temperature 
T unless ctjj = 1 and fyj = ±1 for all pairs ij. 

From Eqs. (3T6)-(3T9) and the expressions of Table I one gets the following 
estimates: 

n~ 1 J i j[m i Vi} = 0, n^ 1 JijlhUi) = -fij^ai (<Jifli + Ujfij) , (5-4) 



-Pi 



"TOt 



-rot 



(5-5) 



rot 



'J 



rriiKjTl 



ft, 



rot 



27™ + -miKmajfli ■ Hj - fa 



+-o-iO-jf2i • flj 




(5-6) 



where we have introduced the effective collision frequency 



4V27T 



-XijrijVij' 



+ 



3 

mi 



(5-7) 



Equations (5-5) and (5-6) are the main results of this paper. They express the 
collisional rates of change as functions of the local values of raj, rij, fii, f2j, T* r , Tf ot , 
T- r , and TJ ot , as well as of the mechanical parameters raij, rrij, <Ji, <jj, Kj, kj, aij, 
and fy. The energy production rates (5-5) and (5-6) can be decomposed into two 
classes of terms. The first class is made of terms headed by 1 + and 1 + fyj which 
exist even if collisions are conservative. These terms do not have a definite sign, are 
proportional to temperature differences (except for the scalar product Q; L • f2j), and 
tend to make the four temperatures equal . The second class is made of terms headed 
by 1 — afj and 1 — j3f- and are positive definite, thus contributing to a decrease of the 
temperatures due to energy dissipation. The terms of the first and second classes 
can be termed equipartition and cooling rates, respectively. Only the latter class 
contributes to the net cooling rate defined by Eq. (5-3). The result is 



4nT 



m. 



n 



(1 




0; 



T~TOt 

i 



+ 



T~TOt 

3 



+ -o-io-jfii ■ rij 



T 



tr 



+ 



T->tr 
3 



(5-8) 



In the case of smooth spheres 
obtained in Ref. 7). 



T), Eqs. (5-5) and (5-8) reduce to those 
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Before closing this section, note that in the monodisperse case Eqs. (5-5), (5-6), 
and (5-8) become 



* 2 



2 K i, o2\ « m 2 A T rot + Kma 2 ^ 2 /12\ 



rot _ " 1 + f> T 



tr 



-rot 



21 + k T rot 

T rot + w 2 fl 2 /12 



x 1- 



c 



£tr ytr _|_ ^rot^rot ^ jitr 

ytr _|_ ^rot 2 7^ tr -|- yrot 



+ Kmo- 2 I? 2 /12 k 



a-/» t^^-tt^ (1+/3) 



(5-9) 



(5-10) 



2 1 - /3 2 / T rot + Km<T 2 ^ 2 /12 
1 — a H — /c + 



1 + k V Ttr 

(5-H) 

where i/ = (4v / 27r/3)x^c r2 v / 2^ tr / ?Tt - Equations (5-9)-(5-ll) agree with those previ- 
ously derived in Refs. 30), 33) with 17 = 0. 

§6. Application to the homogeneous cooling state 

In the so-called homogeneous cooling state (HCS) the flux term V-nj(vj^>(vj, bJi)) 
in Eq. (3-10) is absent. Therefore, the evolution equations for the total and partial 
temperatures are 

rptr rpir rpxot rprot 

d t T = -(T, 5t ^ = -(4 tr -C)^, ft-jr = - (4 rot - C) -jr. (6-1) 

After a certain transient time, a scaling regime is reached where all the time depen- 
dence occurs through the total temperature T, which implies constant temperature 
ratios and equal production rates, i.e., 

£? = & = -=&* = $* = -■ (6-2) 
Since the HCS is an isotropic state, it follows that, for symmetry, i?j = 0. 
6.1. Monodisperse system 

In the monodisperse case the HCS condition £ tr = £ rot yields the quadratic 
equation 

1 — K K frprot rptr \ 



1-a 2 -- — -(1-/3 2 )-— — — (1 + /3) 2 — =0. (6-3) 

1 + k 1 (i + k) 2V p; \T tr r ro V v 7 

In the smooth-sphere limit (ft — > — 1) we get 

T rot J , il±Sii( 1 _ a 2)( 1+/3 )-2_ >00) a<1> 



Ttr ^l^u+z 3 )^ « = i- (6 ' 4) 

Thus the elastic-sphere limit (a — > 1) and the smooth-sphere limit (ft — > — 1) do 
not commute. If the spheres are inelastic (a < 1) but perfectly smooth (ft = — 1), 
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Fig. 2. HCS of a dilute equimolar binary mixture with an = an = Q22 = a, fi\\ = P12 = P22 = /3, 
ki = K2 = |, (72/ 'ci = 2, and mi/mi = 8. Plot of the temperature ratios (a) T* 1 /T{ ot , (b) 
T 2 rot /Ti ro \ and (c) T? /T? vs 1 + /3 for [from top to bottom in (a) and from bottom to top in 
(b) and (c)] a = 1, 0.9, 0.8, 0.7, 0.6, and 0.5. The circles at 1 + j3 = in (c) represent the 
results obtained in the case of perfectly smooth spheres (rotational degrees of freedom ignored 
from the beginning). 7 ' 



the rotational and translational degrees of freedom are decoupled and T rot does not 
change with time, while T tr keeps decreasing due to inelasticity. 34 ) As a conse- 
quence, the ratio T TOt /T tT diverges in the long-time limit. On the other hand, if first 
we assume that the spheres are perfectly elastic (a = 1) and then consider small 
roughness (f3 — > —1), the coupling between T rot and T tr is weak but not broken. 
As long as T tr ~ ^rot^ translational temperature decays more slowly than the 
rotational temperature (£ tr /£ rot ~ k), resulting eventually in a temperature ratio 

T rot/ T tr ^ 1 + p _^ q_ 

6.2. Binary mixture 

In the particular case of a binary mixture, there exist three relevant tempera- 
ture ratios that can be chosen in different ways. Here we take one of the transla- 
tional/rotational ratios (T| r /T[ ot ) and the two component /component ratios (Tg /T-f r 
and T2 0t /T[ ot ). The dimensionless parameter space is twelve-dimensional: the three 
coefficients of normal restitution (an, «i2, 022)5 tire three coefficients of tangential 
restitution (/3n, /?i2, fe); the two parameters K\ and K2, the mass ratio mi/m2, the 
size ratio ai/a2, the mole fraction xi = ni/(n\ +17,2), and the total packing fraction 
4> = (7r/6)(nicjf + n 2 crf ). 

To illustrate the results, here we first assume an equimolar mixture where all 
the spheres are uniformly solid and are made of the same material, the size of the 
spheres of one component being twice that of the other component. More specifically, 
xi = 5, a n = a 12 = a 2 2 = a, f3 u = /3 12 = P22 = P, «i = «2 = f , 02M = 2, and 
7712/7711 = 8. Moreover, we consider a dilute granular gas (<j) <C 1), so that Xij ~ 1- 
Thus only the parameters a and /3 remain free. Figure 2 shows the three independent 
temperature ratios as functions of the roughness parameter 1 + /3 for several charac- 
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Fig. 3. HCS of a dilute equimolar binary mixture with an = oiyz = Q22 = a, /3n = /3i2 = ^22 = /?, 
ki — > 0, K2 = |, T2/T1 — 1, and m,2/mi — 1. Plot of the temperature ratios (a) T^ r / (T[ ot / k\) , 
(b) ^/(Tf^/Ki), and (c) T% T /T\*) vs 1 + /3 for [from top to bottom in (a) and from bottom 
to top in (b) and (c)] a = 1, 0.9, 0.8, 0.7, 0.6, and 0.5. 



teristic values of the inelasticity parameter a. As happened in the monodisperse case, 
the translational/rotational temperature ratio Tl r /T[ ot exhibits a peculiar behavior 
in the smooth-sphere limit 1 + j3 — > 0: it diverges for elastic particles (a = 1) while 
it vanishes for inelastic particles (a < 1). This phenomenon has a reflection in the 
rotational/rotational ratio: either T| ot /T{ ot converges to a finite value or it diverges, 
depending on whether a = 1 or a < 1 , respectively. Quite interestingly, the huge dis- 
parity between the rotational and translational temperatures in the smooth-sphere 
limit has a non- negligible effect on the translational/translational ratio T^/T^ 1 " if 
a < 1: it tends to a finite value different from (in fact higher than) the value directly 
obtained in the case of perfectly smooth spheres. 7 ) Thus, a tiny amount of roughness 
has dramatic effects on the temperature ratio T^/T^, producing an enhancement 
of non-equipartition. 

As a second example, we now consider a dilute equimolar mixture of spheres 
externally identical (same mass, size, and coefficients of restitution), except that 
the mass of the spheres in one of the components is practically concentrated in the 
centre (k± — > 0), while the other component is made of hollow spheres («2 = §)• 
Were the spheres strictly smooth, then the system would be indistinguishable from 
a monodisperse system. However, again the smooth case is singular and the results 
show that both components have different temperatures even in the limit 1 + ft — )■ 0, 
as illustrated by Fig. 3. Note that in the limit K\ — > the rotational temperature 
T[ ot vanishes but the ratio T[ ot /ni is well defined. 

6.3. Locus of equipartition 

As is well known, non-equipartition prevails in granular mixtures. This has 
been illustrated in Figs. 2 and 3. On the other hand, by fine-tuning the mechanical 
parameters and the composition, it is in principle possible to reach equipartition, 
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Fig. 4. HCS of a dilute binary mixture with an = a?i2 = «22 = a, Pu = /3i2 = /?22 = P, and 
ki = K2 = «• (a) Loci of values of a and /3 for « = |, |, and where equipartition holds, 
provided the second equality of Eq. (6-6) is satisfied, (b) The hatched region represents the 
values of mi/m2 and and 01/02 where equipartition is possible. The boundaries of the region 
correspond to the extreme limits 711/712 — » and n\/ni — > 00, while the dotted line corresponds 
to 711/712 = 1. 



i.e., Tf T = T/ ot = T. In this equipartition case, Eqs. (5-5) and (5-6) become 



1 j 



m,i + rrij 



q?. + Kij 
l] 1 + K. 



■1 J 



trot 



i-A 



■1 j 



where v. 



4V27T 



y/T(mi + mjj/rhimj. 



(6-5) 

To show under which conditions 



'v ~ 3 

equipartition is possible, let us consider again a dilute binary mixture with an 
ai2 = «22 = Oi, pn = j3\2 = @22 = P, and m = K2 = K - Otherwise, a, (3, n, mi/1712, 
o"i/o"2, and n^jnx are arbitrary. After simple algebra, the HCS conditions (6-2) yield 



a 



1 — K 
1 + K 



(1-/3 2 ), 



n 2 



at, 



7)71+7772 

2m 2 



(7 



12 



7771+7772 

2mi 



(6-6) 



The first equality establishes a relationship between both coefficients of restitution 
and k that is independent of composition, masses, and sizes of the particles [see Fig. 
4(a)]. Moreover, the composition is constrained by the second equality. Without 
loss of generality, let us take rrt\ < 771,2. Then, positivity of the right-hand side of the 
second equality of Eq. (6-6) implies that max {0,^12} < 01/02 < < 1, where 

1/4 



A 
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mi{mi + mj)/2m- 



1. The above inequality defines a wing-shaped 

region in the plane <ti/<72 vs m\/m2 [see Fig. 4(b)] where equipartition is possible, 
provided both equalities in Eq. (6-6) are satisfied. 
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§7. Concluding remarks 

Despite the success and prevalence of the inelastic smooth-sphere model to char- 
acterize granular matter under rapid flow conditions, ^ it is known that the presence 
of friction in collisions can be relevant and even essential for a more realistic model- 
ing of granular flows. 6 )' 41 ) Moreover, granular matter is typically present in nature 
in polydisperse form. In this paper we have combined both aspects (roughness 
plus polydispersity) in a model of inelastic hard spheres with constant coefficients 
of normal and tangential (ftij) restitution. Otherwise, the masses, diameters, 
moments of inertia, coefficients of restitution, and composition are arbitrary. Given 
the complexity of the system, here we focus on the derivation of the (partial) energy 
production rates and ^[° t due to binary collisions. These quantities determine the 
collisional rates of change of the translational (Tj tr ) and rotational (JJ ) granular 
temperatures associated with each component, as well as the cooling rate £ of the to- 
tal temperature of the mixture. They can be considered as the most basic quantities 
in a granular mixture since they are involved in the energy balance equations. 

Starting from the collisional rules (2-6) and (2-7), together with Eq. (2-12), the 
energy production rates are first expressed in a formally exact way in terms of the 
two-body distribution function fj9 [see Eqs. (3-11), (A-7), and (A-8)]. Next, the true 

function is replaced by its pre-collisional orientational average f\- defined by Eq. 
(3-13), what is justified if the density and/or the inhomogeneities are not too large. 
This replacement allows us to express the collisional rates of change as combinations 
of two-body averages, Eqs. (3-16)-(3-19). To get more explicit results as functions of 
the local values of densities, temperatures, and mean angular velocities, a maximum- 
entropy approach is followed in Eq. (4-2). The final expressions are given by Eqs. 
(5-5), (5-6), and (5-8). These are the main results of the paper and generalize previous 
results derived for monodisperse rough-sphere gases 37 ) and polydisperse smooth- 
sphere gases. 7 ) 

As a preliminary application of our results, we have studied the HCS, where 
collisions are the only source of energy change. For a binary mixture of common 
coefficients of restitution = a and ftij = ft, we have analyzed the influence of both 
inelasticity and roughness on the three independent temperature ratios X^" /T{ , 
T2 /T\ ■> an d T2 0t /T[ ot . As a surprising result, we find that a small amount of 
roughness has a significant effect on the translational/translational ratio /T{ T . 
More specifically, it turns out that, at a given value of a, the value of l^i lTl t ne 
limit ft — > — 1 differs from the value obtained directly in the smooth case (ft = —1). 
This paradoxical phenomenon is mainly due to the fact that, even if ft is close to 
— 1, there exists a transfer of energy from the rotational to the translational degrees 
of freedom that becomes relevant since the rotational temperatures are much larger 
than the translational ones in the nearly smooth case. The latter effect is closely tied 
to the non-stationary character of the HCS. In homogeneous steady states, as the 
one driven by a white-noise thermostat, the extreme sensitivity of the ratio T^/T^ 
to whether ft — > — 1 or ft = — 1 disappears. 

The results obtained here can be applied to several problems. First, we plan to 
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compare the theoretical results derived in this paper with computer simulations in 
homogeneous states. This will allow us to assess the reliability of the approximations 
made here to get explicit results. Secondly, taking the local version of the HCS as 
the reference state, a Chapman-Enskog method can be followed to get the Navier- 
Stokes constitutive equations for a mixture of rough spheres. Finally, following steps 
similar to those done in Ref. 46), a kinetic model for rough spheres preserving the 
energy production rates obtained here will be proposed. 
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Appendix 

Collisional changes 

In this Appendix the expressions for the collisional changes of Vj, Wj, v 2 , and to 2 
are given. First, note that from Eq. (2-12) it follows that 

Qij = 5?,K- • a) 2 + [(a x Vij ) 2 + (?x S,,) 2 - 2v t] ■ (a x S tJ )] , (Al) 

a x Qij = fa [a x - ?x(5x S^)] , (A-2) 

(a x Qyf = /3 2 [(a x v^-) 2 + (? x S,,) 2 + 2(5 x v tj ) ■ S -] . (A-3) 
where use has been made of the mathematical identities 

a x (<r x A) = (a ■ A)a - A, (5 x A) • (5 x B) = A • B - (5 • A)(a ■ B). (A-4) 

Next, Eqs. (2-6) and (2-7) yield 

rriiv'i - rriiVi = -5ij(v^ • a)a - [v^- - ax - (v^ • a)a] , (A-5) 

- IiUi = ~Aj [S- x Vij ax (a x S^)] , (A-6) 



Mixtures of inelastic rough hard spheres 



17 



/2 2 

m,iV i — rriiVi 



2 - B 2 

^(Vij • £) 2 + [(o 1 x Vij ) 2 + (a x S -) 2 - 2v - • (5 x Sij)] 



mi nii 
-2aij(vij ■ <r)(v, ■ a 



28ij [(a x Vj) • (<r x Vy) -v r (?x S^)] , 

(A-7) 



/ iW J 2 - = — ^" [(£ X Vij ) 2 + (5 X Stf) 2 + 2(5 X Vy) • Sy] 

TTliKi 

-PijUiUi ■ [a x Vy - x (ct x Sjj)] . 



(A-8) 



Similar expressions are obtained for particle j by exchanging i f> j and <r <h> — er. 
Combining Eqs. (A-7) and (A-8), plus their counterparts for particle j, one can 
get Eq. (2-14), where use is made of Eq. (2-11) and the identity A • (<r x B) = 
-B • (<r x A). From Eqs. (A-5)~(A-8) one can easily obtain Eqs. (3-16)-(3-19) by 
using the mathematical identities 



/ 



2tt , 

v- ■ 

£+1 v 



/ 



J da0(vij ■ ff)(vjj • a) 1 a a 



2tt 



(£ + !)(£ + 



(A-9) 
(A-10) 
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where I is the unit tensor. 
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